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ABSTRACT 


Perfect liquid mixing but no solid diffusion is assumed in 
normal freezing. In addition, the molar compositions of the freez- 
ing solid and remaining liquid, respectively, follow the solidus 
and liquidus curves of the constitutional diagram. For the linear 
case, in which both the liquidus and solidus are perfectly straight 
lines, the normal freezing equation giving the fraction solidified 
at each melt temperature and the solute concentration profile in 
the frozen solid was determined as early as 1902, and has since 
been repeatedly published. Corresponding equations for quadratic, 
cubic or higher-degree liquidus and solidus lines have also been 
obtained. This memorandum gives the equation of normal freezing 
for ideal ternary liquid solutions solidified into ideal solid 
solutions of the pseudobinary type. Sample computations with the 
use of this new equation were made and are given for the Ga-Al-As 
system. 
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INTRODUCTION 


Perfect liquid mixing but no solid diffusion is assumed in 
normal freezing. In addition, the molar compositions of the 
freezing solid and remaining liquid, respectively, follow the 
solidus and liquidus curves of the constitutional diagram. 

For the linear case, in which both the liquidus and solidus 
are perfectly straight lines, the normal freezing equation giving 
the fraction solidified at each melt temperature and the solute 
concentration profile in the frozen solid was determined as early 
as 1902 (Ref. 1), and has since been repeatedly published (Refs. 2-9). 
Corresponding equations for quadratic, cubic, or higher-degree 
liquidus and solidus lines have also been obtained (Refs. 10-12). 

A recent publication presents the normal freezing equation for 
ideal dilute solutions or alloys in which the liquidus and solidus 
are exponential functions of the reciprocal absolute melt tempera- 
ture (Ref. 13). The logarithm of the mole fraction segregation 
coefficient for the dilute component is, therefore, also a 

linear function of the reciprocal of the absolute temperature. 

Further, in such ideal solutions, the ratio of partial pressure 
of each component, p^, to that of the pure component, p°, 
equals the atomic fraction x^, i.e., x^ = p^/p°. This last 

equation permits the computation of the solidification behavior 
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of such solutions or alloys and the resultant solid concentration 
profiles from basic thermodynamic quantities without the use of 
phase diagrams . 
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THE PSEUDOBINARY SYSTEM 


This memorandum gives the equation of normal freezing for 
ideal ternary liquid solutions of A, B, and C atoms (Fig. 1) 
solidified into ideal solid solutions of AB-AC compounds with 
A atoms on one sublattice and B or C atoms on the other. 
Such systems can be treated as pseudobinary systems. Account is 
taken of the entropy of mixing of the three atomic species to 
form the ideal liquid solutions along the pseudobinary. 



Fig. 1 Constitutional Diagram for a Ternary 
System of the Pseudobinary Type 
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( 1 ) 


For the compound AB, we have 

AB(s) = A(l) + B(l) . 

The equilibrium constant for the reaction is 



( 2 ) 


and the standard free energy AF° B and the free energy of forma- 
tion AF^ B are related to K at temperature T by 

■Kb + - rt ln k ’ < 3 > 

For the compound AB, the free energy of formation, AF^ B> 

is related to the free energy of melting, AF™ , and the heat of 

AJj 

melting, at the melting point by 

af ab - ah ab (1 ■ t^ab 5 - -Kb - RT ln 4 • <4) 

It has been assumed that the difference in heat capacity be- 
tween the liquid and solid is negligible. The -RT • ln 4 term 
is the ideal entropy of mixing to form the stoichiometric liquid 
phase. Hence, 

1 # 1 

-RT ln 4 • - AH^ b (1 - T/t” b ) . (5) 

X AB 

Similarly, for the compound AC, we have 
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1 . 1 

- RT ln 4 • "*T-k ' ah ac (1 ’ T/T Ic ) > 

X AC 


( 6 ) 


where AH^, is the heat of formation of the compound AC, assumed 


„m 


constant; and T. r is the melting point of AC. 

For the ternary liquid phase of the pseudobinary, 


2 » 


1 ^ 1 1 
Xb + X C - J 


(7) 


Hence , 


-RT ln 4 • b 


and 


- AH AB <! - T/T Ib> 


( 8 ) 


■RT ln 4 • i ■ ( i - x^)/(X - x® B ) 


ah ac (1 - t/t ac> ' < 9 > 


Equations (8) and (9) are essentially the same as those previously 
given (Refs. 14,15). From these equations, 


2 ' *B ■ X AB ' exp 


AH 


m 


AS N /t . i/r m 
R l 1/T 1/T AC 


\ 


= 4 b • exp(a + b/T) = x= B • U 


(10) 


1 ' 2 • x B = (1 - x® B ) 


= (1 ’ X AB> 


ah ac 

[i/T - 1/T“ C ] 

R 

1 

j 


exp 


exp(f + g/T) = (1 - x Afi ) 


(ID 


V 


where 
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( 12 ) 


b '- AH AB /R 1 3 " b/T AB ' S " AH AC /R ’ and f = ‘ g/T AC 

are constants for a given ternary system of the type shown in 
Fig. 1, while the temperature-dependent variables 

U = exp(a + b/T) and V = exp(f + g/T) . (13) 

Equations (10) and (11) yield the solidus equation for the 
compound AB 

1 - x !b ' U ■ (1 • X AB> ' V • <14) 

Hence , 

X AB = (1 ’ V)/(U ‘ V) • (15) 

g 

The value of x^ can be as small as zero when, according to 
Eqs . (15) and (13), V = 1 = exp(f + g/T), i.e., when T = -g/f = 
T™ c and x^ = x^/2 = ® [Eq. (H) 3 » as expected. The value of 

g 

x.., can also be as large as 1 when U = 1 = exp(a + b/T) in 
Ad 

accordance with Eqs. (15) and (13), i.e., when T = -b/a = T™ fi 
1 s 

and x B = x AT1 / 2 = \ [Eq. (10)], also as expected. 
d AB 

Notice that there are no AB or AC compounds in the melt, 
while in the solid there are only compounds AB and/or AC, but 

no individual elements A, B, and C. Hence, the thermodynamic 

, ^ . s s 1 1 1 

molar concentrations X ^ B > X AC’ X A’ X B’ X C ' 

Considering the ternary liquid or solid solutions as a whole, 

and, without regard to crystallographic lattice arrangements, we 
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can define the chemical molar concentrations 


V V 


and m . 


Evidently, 


m A = m A = 2 


1 1 s , s i 

m B + “C = “B + “C = * 




1 

“b 




and 


m 


B 


= * X AB 


( 16 ) 
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THE NORMAL FREEZING EQUATION 


We can now balance the amount of component B in the remain- 

1 s 

ing liquid by its chemical molar concentrations m^ and m^ 

then obtain the usual differential equation for normal freezing 
(Refs. 9, 11-13) as follows 


dm B 


1 s 
"b ' m B 



(17) 


where p is the fraction solidified. 

Substituting into the normal freezing differential equation 
results in 


where 


and 


- d i”d - p) - I v/ u ) d(1/T > 

. [ a . _ - P _ , __L_J Af ±\ 

1(1 - U) (1 - V) (1 - v/u)j V 



(18) 
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and 


q = exp(f - ga/b) '' 

ag - bf 

r = exp — ° — 

* g “ b 

bf - ag 

s = exp “ — & \ 

U c = exp(f + g/T o ) 

V Q = exp (a + b/T Q ) 


( 19 ) 


_ 1 

1 I-* 

1 1 

c: 

|o 

a 

i— * 

i 

< 

o 

■ 

(3 

[1 - V /U 
o o 

p = 1 - 

11 - u J 


[1 - V j 


[l - V/U 


( 20 ) 


where T q is the initial freezing temperature and can be deter- 
mined from the pseudoliquidus line on substituting for x^ fi the 
initial melt composition, x^ g q - 2 • x^ q . 
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THE Ga-Al-As SYSTEM 


This ternary system can be considered as a GaAs-AlAs pseudo- 
binary. In the liquid melt, we have 


x . = x. 

A As 


"As " °' 5 


1 , 1 1.1 1.1 1.1 nc 

x_ + x_ = x_ + x A1 ® m_, + m_ = m_ + m Al = 0.5 
B C Ga Al HB C Ga Al 


In the GaAs-AlAs mixed crystal, we have 


s 0 s 

x„ A = 2m_ 
GaAs Ga 


s 0 s 

X AlAs 2 Al 


s s 

X GaAs X AlAs 


2( ”Ga + “Al 5 ’ 1 ‘° 


To determine the phase diagram [i.e., liquidus and solidus, 
Eqs . (8) — (15) ] and, in addition, to compute the normal freezing 
behavior according to Eqs. (18) -(20), the following four thermo- 
dynamic constants are required (Ref. 14 gives the values of these 
constants) : 


Melting point of GaAs 
Melting point of AlAs 
Heat of formation of GaAs 
Heat of formation of AlAs 


1511 K 
2013 K 

16.64 Kcal/mol 
22.8 Kcal/mol 
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Hence, a = 5.700, b = -11,470, f = 5.542, g = -8,374, q = 3.983, 
r = 0.1505, s = 166.5, a = -1.335, £ = 0.1772, and y = 1.006. 

Table 1 summarizes the solidification behavior of the Ga-Al-As 
system, treated as a pseudobinary. The first three columns describe 

the phase diagram, i.e., the concentration of aluminum in the liquid 

1 s 

melt and that of AlAs in the solid crystal at tem " 

peratures from 1720°C down to 1238°C at 20°C intervals. The 

next two columns, obtained from the first three, show that the 

s 1 

segregation coefficient k for Al, i.e., = ^ias'^AI* 

g 

always exceeds unity while that of Ga, i.e., k^ a = (1 - *- AlAc )/ 


AlAs' 


(1 - 2x a ^) , never does so except at 1238°C. Both of these segre- 
gation coefficients, however, are widely variable, increasing with 
decreasing melt temperature, by factors of over 6.3 in the case 

of k.. and of over 3.7 in the case of k_ . Thus, for this 
Al Ga 

system, as well as for many others, the linear normal freezing equa- 
tions (Refs. 1-9) involving constant segregation coefficients are 
generally not applicable. 

These phase diagram data can, of course, be obtained also by 
means of the equations of Ilegems and Pearson (Ref. 14), or those 
of Steininger (Ref. 15). The direct computation of the normal 
freezing data from basic thermodynamic constants, however, appears 
to have never been attempted or published. By means of Eqs. (18) 
and (19) of this memorandum, these normal freezing data are cal- 
culated and given as the last nine columns. Here, the fraction 
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SOLIDIFICATION OF Ga-Al-As SYSTEMS 
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Basic thermodynamic constants used in the computation: 

"Ws ■ 16 ' 64 lteal/ " < ’ 1 T C«A 8 ■ 1511 K 

4H A1As ' 22 ' 8 Kcal/ ” 01 T A!As ‘ 2013 K 



solidified, p, is given as a function of both the melt tempera- 


ture, T, and the Al concentration in the original liquid melt, 
q . The concentration profile can be obtained by plotting the 

g 

concentration of AlAs in the solid crystal, x aias* against the 


fraction solidified, p, for the appropriate initial Al concen- 
tration in the melt, x}, 

5 Al,o 
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OTHER SYSTEMS 


Many other ternary systems can also be treated as pseudo- 
binaries. The necessary thermodynamic constants are usually 
available in the literature. Steininger (Ref. 15), for example, 
has examined a number of important systems and given the neces- 
sary thermodynamic constants required to compute the solidifica- 
tion behaviors and the phase diagrams. These systems include: 
InAs-GaAs, InSb-InAs, InSb-GaSb, GaSb-AlSb, InSb-AlSb, CdTe-ZnTe, 
CdTe-CdSe , HgTe-CdTe, GeTe-MnTe, PbBr 2 -PbCl 2 , InAs-InP, SnTe-PbTe, 
PbTe-PbSe, HgTe-HgSe, ZnTe-ZnSe, and HgTe-ZnTe. For the first 
13 systems, his calculated solidus curves are in good to excellent 
agreement with the published experimental values. These systems 
can, therefore, be considered ideal. One can thus compute with 
some confidence their normal freezing behaviors according to 
Eqs . (18) -(20) of this memorandum. 
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